The electron spectrum of epitaxial graphene monolayers 
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Epitaxial graphene on SiC possesses, quite remarkably, an electron spectrum similar to that of 
freestanding samples. Yet, the coupling to the substrate, albeit small, affects the quasiparticle 
properties. Combining ab initio calculations with symmetry analysis, we derive a modified Dirac- 
Weyl Hamiltonian for graphene epilayers. While for the epilayer on the C-face the Dirac cone 
remains almost intact, for epilayers on the Si-face the band splitting is about 30meV. At certain 
energies, the Dirac bands are significantly distorted by the resonant interaction with interface states, 
which should lead to mobility suppression, especially on the Si-face. 
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Continuing improvement of epitaxial graphene on 
SiC that culminated in observations of the quantum 
Hall effect 5,4 and of a very high mobility at the Dirac 
point have raised the ranking of this material dramati- 
cally. Apart from the capability to emulate freestanding 
graphene, epitaxial graphene has a number of specific 
qualities making it an intersting material in its own right. 
These features stem from the interaction with the under- 
lying substrate - the Si or C terminated SiC surface. In- 
terestingly the growth of graphene differs drastically on 
both surfaces. On the Si-face the growth is slow thus facil- 
itating the fabrication of a single monolayer. 1 The subse- 
quent C-layers arrange in the common graphite-type AB 
(Bernal) stacking. On the C-face, the much faster growth 
typically yields multilayer stacks of mutually rotated C- 
layers. 2 The rotation is very important since it decou- 
ples individual layers electronically such that the whole 
stack behaves effectively as a single graphene sheet. 2,6 
Recently also monolayer graphene has been achieved on 
the C-face. Compared with the C-face, the electron mo- 
bility in Si-face graphene is much lower. Along with the 
preference for Bernal stacking of the latter, this indicates 
a stronger coupling to the substrate. 

In fact, the graphene monolayer does not reside 
directly on SiC but rather on some buffer layer, 
as first realized theoretically 8,9 and then confirmed 
experimentally. ' The currently accepted buffer model 
for the Si-face is a corrugated graphene layer, that 
is covalently bonded to the substrate fitting into a 
(6-\/3x 6\/3)R.30 surface reconstruction. 1 1 The reconstruc- 
tion unit cell almost exactly coincides with a (13x13) 
graphene unit providing the commensurable base for sub- 
sequent graphene layers. The strong covalent interaction 
with the substrate completely erases all Dirac- Weyl fea- 
tures of the buffer. Hence it is the second carbon layer, 
which exhibits the graphene-like band structure and is re- 
ferred to as "monolayer graphene" . 8,!1 This scenario was 
convincingly confirmed in recent experiments with hydro- 
gen intercalation. 1 " Diffusing underneath the buffer H 
atoms cause its release and as a result, the formation of 
a quasi-freestanding Bcrnal-stacked bilayer was observed. 

The situation on the C-face is less clear. A distinct 



buffer layer has not yet been identified and it was specu- 
lated, that already the very first carbon sheet might be 
graphene-like. However, this contradicts calculations 8,11 
which clearly show the extinction of the Dirac spectrum 
of the first C-layer on both SiC surfaces. Si adatoms 12 or 
a corrugated C-layer 8 ' 13 were considered as buffer mod- 
els on the C-face. Generally, the much weaker graphene- 
substrate coupling on a C-face suggests, that in this case 
the particular interface structure is not as important as 
on the Si-face. 

The graphene-substrate coupling, and especially its ef- 
fect on the electron spectrum in the vicinity of the Dirac 
point, is of particular interest for electron transport. It 
has been a subject of a long debate whether the graphene- 
substrate interaction opens an energy gap with experi- 
mental estimates ranging from zero to 0.3 eV. 1,14,15 

In this article we consider the electronic structure of 
graphene monolayers on SiC combining ab initio calcula- 
tions with symmetry analysis. We derive the low energy 
Hamiltonian, that replaces the Dirac- Weyl Hamiltonian 
of a freestanding graphene. Instead of the computation- 
ally demanding (6-\/3 x 6-\/3) structure, we adopt a strain 
free, commensurable (5x5) interface model with a corru- 
gated carbon layer as a buffer (cf. Fig. 1). The (5x5) 
reconstruction is observed, although more rarely, on the 
Si- face. 111 With the same buffer model for the C-face we 
find a much weaker graphene-SiC coupling than on the 
Si- face. 

Backfolding of the graphene K and K' points to the 
T-point in the (5 x 5) Brillouin zone produces four closely 
lying energy branches. On the Si-face, we find that the 
Dirac cone is split by about 30meV, whereas on the C- 
face two branches of an essentially unperturbed cone exist 
and the other two are separated by a very small gap 
(< lOmeV). The splitting of the Dirac cones is not due 
to the corrugation of the graphene layer, but arises from 
the interaction with the substrate. We verified this by 
calculating the freestanding carbon layer with the same 
atomic positions, for which we found a perfect Dirac cone. 

The buffer layer is commensurate with the SiC surface 
when rotated by ~ ±16.1 °, such that the resulting struc- 
ture possesses a (5x5) periodicity relative to the SiC sub- 
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strate (cf. Fig. 1). The following graphene epilayer can 
be either rotated by ~ =f32.2 ° with respect to the buffer 
or aligned in AB-stacking. As seen in the right panel 
of Fig. 1, the (5x5) structure naturally results from a 
commensuration of the twisted carbon bilayer and the 
SiC surface. The (5x5) SiC unit cell is almost perfectly 
commensurate with the (\/l3X\/l3) graphene cell. No- 
tably, the rotation angle of (30 ± 2.2) ° dominates in the 
multilayered graphene stacks 2 and a 30 °-rotation of the 
graphene bilayer relative to the SiC cell is a common fea- 
ture for both Si- and C-terminated surfaces. > 




FIG. 1. (Color online) Top and side views of the rotated 
graphene monolayer on a (5x5) SiC(0001)/buffer slab with 
the H-passivated bottom (C/Si atoms: small dark/large 
light spheres). The shaded area indicates the (V13x\/T3)- 
graphene Wigner-Seitz cell of the graphene/buffer bilayer 
(medium-sized/small spheres) with relative rotation of 32.2 °, 
that is commensurate with the (5x5) SiC cell (surface atoms: 
large spheres) - the three central hexagons of the graphene 
layer are emphasized along with the vector t (cf. text). Lower 
left: The Brillouin zone of an unperturbed graphene layer. 

We used the VASP density functional package 17 to 
obtain the interface geometry by relaxing atomic posi- 
tions in the (5x5) unit cell and to calculate the band 
structures. The buffer layer shows a significant corruga- 
tion (cf. Fig. 1), that is partly transmitted into the top 
C-layer. The calculated large-scale band structure is de- 
picted in Fig. 2(a)-(d). The graphene-type linear bands 
appear with the second carbon layer [Fig. 2(c), (d)]. The 
buffer does not possess that feature, but supplies a num- 
ber of flat interface states, which pin the Fermi level Ep 
[Fig. 2(a), (b)]. Owing to Fermi level pinning, the epilayer 
is cither n- or p-doped on the Si- and C-face respectively, 
in accord with the measurements. 1,2,19 It is also visible, 
that the resonant interaction with the interface states 
causes a significant distortion of the Dirac bands. 

Figures 2(e)-(j) show the high-resolution energy spec- 
tra close to the Dirac point. Regardless whether the 
twisted or the Bernal-type structure is chosen, the en- 
ergy spectrum resembles a modified single layer spec- 
trum (with folded K and K' points) and not a graphene 
bilayer spectrum. This again confirms, that the buffer 



does not possess the Dirac states at the K-point, i.e. it 
is completely passivated by the substrate. The substrate 
potential is rather weak, especially on the C-face, caus- 
ing the band splittings in the range of 5 — 30 meV. As 
pointed out above, the slight warping of the epilayer does 
not affect its band structure within an accuracy of 2 meV. 
Hence it should be possible to obtain the low-energy spec- 
tra by accounting for the substrate as a perturbation of 
the ideal graphene states. In the following we use symme- 
try considerations to construct the effective Hamiltonian 
for the graphene monolayer. The derived Hamiltonian re- 
produces the calculated ab initio energy spectra and can 
serve for the description of the quasiparticle dynamics. 

Commonly, the graphene spectrum is introduced via 
the tight-binding modeling of the 7r-bands. 2 " Yet, the 
symmetry underlying the spectrum is much more trans- 
parently expressed within the empty lattice approach. 
The three states |Kj) (i = 1,2,3), that correspond to 
the plane waves <px i (r) at equivalent corners Ki of the 
Brillouin zone, combine into the K-point eigenstates and 
similarly for the time- reversed states at points K[. The 
trigonal symmetry dictates the Hamiltonian matrix 
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where the mixing of the basis states is expressed via the 
symmetrized operator of a ( 27r /3)-rotation: 83 + Cg = 
exp (i^-L z ) + h.c, with the angular momentum opera- 
tor L z in the Hilbert space 1 = 1. The diagonal elements 
in Eq. (1) are linear invariants of a small displacement 
k from the K-point and the unit vectors rii — K,-/|Kj|. 
The phase ip of the matrix elements V = |V| exp (itp) de- 
pends on the choice of the coordinate origin. For ideal 
graphene (f takes the values or ± 27r /3. In an epitax- 
ial layer, however, the inversion symmetry and hence the 
mirror planes of the small group C% v are lost. This im- 
plies a reduction to C3 and allows arbitrary values of (p. 
By time reversal the Hamiltonian at the iC-point is 
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Replacing L z by its eigenvalues (0,±1) we obtain the 
energy levels at the K-point 

e = 2\V\cosip and e± = 2\V\ cos f~ ± cpj . (3) 

Of special interest are the two states, that form the tip of 
the Dirac cone in ideal graphene. Depending on ip these 
can be any two of the levels in Eq. (3). The choice ip = 
for the coordinate origin in the center of the graphene 
hexagon selects e+ and £_. 

The states |Ki) and |K^) are coupled by the Umklapp 
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FIG. 2. (Color online) Calculated band structure of epitaxial graphene on SiC(OOOl) (C-face, upper row) and on SiC(OOOl) 
(Si-face, lower row). Large-scale energy spectra (a),(b) of the buffer on SiC and (c),(d) of the 32.2° rotated buffer/epilayer 
structure. The high resolution spectra in the range p/20 KT, 1 /20 TM] are shown for (e)-(h): the rotated buffer/epilayer with the 
equilibrium buffer-epilayer distance [(e), (f)] and with a distance artificially increased by 1 A [(g), (h)]; Bernal stacking [(i),(j)]. 
The color code indicates the degree of localization of states on the graphene epilayer (pa) and on the buffer layer (pb) as 
deduced from the integrated density g(z) = J \^>\ 2 dxdy. Fitted bands are shown as solid curves (cf. Eq. (9) and Tab. I). The 
insets illustrate the possible types of band structures that follow from Eq. (9). 



process due to the substrate potential V(r). Since we 
account for V(r) perturbatively, we can assume that it 
preserves (within the intervalley matrix elements) point 
symmetry of a pristine surface. However, the point group 
centers of the epilayer and of the substrate (i.e. of the last 
Si or C atomic layer) are displaced by the vector r, which 
connects two adjacent graphene atoms (cf. Fig. 1). This 
displacement generates a phase factor of the intervalley 
matrix elements 

Vij = J ip* Ki (r) V(r + t) <p Ki (r) dr 

= e-< K *- K 'i) r J y* Ki {r)V{r)y K ,{r)dr. (4) 

The shift of the coordinate origin in Eq. (4) enables us 
to exploit the trigonal symmetry, which requires | V13 j = 
\V 2 x\ = \V 32 \. 

Given r and Ki — K'j (cf. Fig. 1) one finds, that the 
phase in Eq. (4) takes values ± 27r /s or 0. For and 
Vji the phase factors can be made complex conjugated 
by shifting phases of the basis K-functions: |Ki) — ► 
e*T [Ki), |K 2 ) |K 2 ), |K 3 ) e-** |K 3 ). This, how- 
ever, changes the phases in Eq. (1), such that the K- 



Hamiltonian takes the form 
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After this phase transformation the intervalley interac- 
tion Eq. (4) acquires a structure similar to Eqs. (2) and 
(5): 



V = a + 2(3cos ( —L z + — 
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In the diagonal representation (L z = diag(— 1, 0, 1)) the 
Hamiltonian matrix reads 
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where p = vp(k x + iky) with the Fermi velocity up- Note, 
that the phase shift in Eq. (5) leads to a cyclic permu- 
tation of the K-point eigenvalues. The "Dirac part" Hjj 
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of Eq. (7) comprises the energy levels, which merge for 
(p. = 0. Introducing parameters A = 2v3|V^| sin ip and 
b = a /2 + (3 we obtain 

/ -A p* 2b \ 

Hd = _ A _ p « ■ (8) 

\ — p A / 

The eigenvalues of ifo are 

= ±|6| (±) ^(ATlfeD' + ^Ffc) 2 (9) 

with i = 1, ...,4 for different combinations of signs. 
Equation (9) describes the four bands originating from 
the K and K' states of an ideal graphene layer subject 
to a symmetry lowering potential of the substrate. The 
spectrum of Eq. (9) takes qualitatively different forms 
for A > \b\, A = |6| and A < |6| (cf. insets Fig. 2). By 
fitting Eq. (9) to the ab initio bands, we determine the 
Hamiltonian parameters as given in Tab. I. The 32 °- 
twisted buffer /epilayer pair on the Si-face and the Bernal 
stacking on the C-face matches the A > \b\ case, whereas 
the twisted C-termination and the Bernal stacking on the 
Si- face correspond to the A = |6| case. 



buffer/ epilayer 


32.2°-twisted 


Bernal- 


-type 


[meV] 


\b\ 


A 


|6| 


A 


Si-face (0001) 


12.7 


21.7 


13.5 


11.0 


C-face (0001) 


2.1 


2.1 


20.3 


29.9 



TABLE I. Parameters A and |i>| as obtained from a fit 
of Eq. (9) to the ab initio band structure. The fit yields 
v F = 0.8«p cc (Si-face) and v F = 0.96 «p oc (C-face)' with 
Vp cc = 8.37 ■ 10 5 ms _1 being the calculated Fermi velocity 
of freestanding graphene. 



The fitted bands of Eq. (9) are shown in Fig. 2(e), (f) 
and (j),(i). A noticeable deviation of the fit from the 
ab initio data occurring for some bands originates in the 
repulsion by the closely lying interface states (see below) . 

Setting A = in Eq. (8), Ho becomes similar to the 
Hamiltonian of the freestanding twisted bilayer, where 



the spectrum is formed by the K-states of the two 
graphene layers. Due to the inequivalence of neighboring 
atoms, the intervalley scattering in Eq. (6) has the same 
structure as in the case of the "odd sublattice exchange" 
of the twisted bilayer. 21 

In graphene epilayers the quasiparticle dynamics is de- 
scribed by the effective Hamiltonian Hp, which replaces 
the Dirac-Weyl Hamiltonian of ideal graphene. An impor- 
tant effect, apparent in Fig. 2, is the strong resonant inter- 
action of the Dirac bands with the interface states. This 
interaction should mediate a quasiparticle scattering by 
interface phonons. By artificially increasing the top layer 
separation we can trace the transition from an epitaxial 
to a freestanding graphene spectrum [cf. Fig. 2(g),(h)]. 
The band splitting as well as the resonant interaction 
decrease sharply, yet, the typical pattern of the resonant 
coupling is clearly seen in Fig. 2(h). Moreover, it is appar- 
ent in Fig. 2(h), that the resonant interaction is subject 
to certain selection rules: while one Dirac branch strongly 
couples with the interface state, the interaction matrix 
element vanishes for the other. Hence, the dispersion 
of only one Dirac band is significantly affected at equi- 
librium separation of the top layer [cf. Fig. 2(f),(i),(j)]. 
The resonant interaction is much more pronounced on 
the Si-face, where the coupling matrix element is about 
50meV [as estimated from Fig. 2(f)]. This is in accord 
with the lower carrier mobility observed on the Si-face in 
comparison to the C-face. 1,2,5 As visible in Fig. 2(c), (d), 
the resonant coupling also occurs close to the Fermi en- 
ergy. 

For a (6v / 3x6 v / 3)R30 SiC(0001) surface a similar 
substrate-mediated interaction of the Dirac bands has to 
be considered, albeit with a different translational sym- 
metry. Due to the weak coupling of the graphene layer 
one expects a qualitatively similar scenario. Indeed, the 
gap opening as well as the formation of flat interface 
"midgap states" were found by ab initio modeling of 
photocmission spectra for a (6-\/3x 6y / 3)R30 substrate. 
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